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Experiments in the antiferromagnetic phase of a quasi 2D 5 = 1/2 quasi-square lattice antifer- 
romagnet Cu(pz)2(C104)2 reveal a biaxial type of the anisotropy, instead of the easy-plane one, 
considered before. The weak in-plane anisotropy, found by means of electron spin resonance spec- 
troscopy and magnetization measurements, is about an order of magnitude weaker, than the off-plane 
anisotropy. The weak in-plane anisotropy results in a spin-flop phase transition for the magnetic 
field aligned along easy axis, and, thereby, in a bicritical point on the phase diagram. A remarkable 
feature of the weak in-plane anisotropy is the abrupt change of its sign at the spin-flop point. This 
anisotropy switching disappears at the tilting of magnetic field to the easy axis by the angle of 10° 
within the plane. The nature of the abrupt anisotropy reversal remains unclear. The phase diagram 
is characterized by the increase of the ordering temperature Tjv in the magnetic field used, except 
for a dip near the bicritical point. 

PACS numbers: 75.40.Gb, 75.50.Ee, 76.50.+g 



I. INTRODUCTION 

Heisenberg 5 = 1/2 antiferromagnet on a square lat- 
tice (HSLAF) is a popular model of low-dimensional 
magnetism^. An ideal HSLAF has no long-range or- 
der except at T = 0, where a Neel-type ground state 
with 40% reduction of ordered spin component should 
be realized 2 -. In real quasi-2D antiferromagnet a weak 
interlayer interaction is present, providing a Neel order 
at T > 0. An organometallic compound Cu(pz) 2 (0104)2 
(copper pyrazine perchlorate) has been considered as an 
example of quasi-2D HSLAF: copper ions carrying spin 
5 = 1/2 arc bridged together in a slightly distorted 
square lattice layers by pyrazine (C4H4N2) rings as shown 
on Fig. [TJ 

Copper pyrazine perchlorate crystallizes from a solu- 
tion in the space group C2/m at room temperature, how- 
ever, at cooling, near 180 K, there is a phase transition 
to a structure which has the space group C2/c i. In t he 
low temperature phase, parameters b and c of the mon- 
oclinic lattice are close to each other. The rectangles 
with sides b and c are approximately squares. Diagonals 
of these rectangles form a rhombic lattice, these rhom- 
buses are slightly distorted squares. The angle between 
the a and c axes differs only a little from 90°, thus the 
lattice may be represented as weakly distorted tetrago- 
nal lattice. Magnetic ions Cu 2+ (5=1/2), placed at the 
corners of rhombuses, form layers in be planes. The near- 
est neighbor exchange paths are symmetrically equivalent 
and exactly identical. Therefore the exchange network 
within the be planes is equivalent to that of a square lat- 
tice. Due to CIO4 complexes between the layers, they are 
well separated, as well as due to half a period in-plane 
shift between layers. Because of this shift a magnetic ion 
within a layer is equidistant from four ions in adjacent 
layer, therefore the interlayer coupling is canceled in the 




FIG. 1. (Color online) Cu(pz) 2 (C10 4 )2 structure. Two lay- 
ers, each containing a square magnetic lattice, are displayed; 
CIO4 complexes are not shown for clarity. Colors of ions and 
connecting lines in the lower layer are faded out. Crystallo- 
graphic data are taken from 0]. 



first order—. 

Indeed, estimation of the interlayer effective ex- 
change Ji from values of Tn = 4.25 K and nearest- 
neighbor exchange J = 18.1 K by an empirical relation 
J±/J ~ cxp (— ^j^-) derived from quantum Monte- Carlo 
simulation^, leads to a very small value of Jj_ — 9 ■ 10 -4 J 
Q. Magnetic moment per Cu 2+ ion in the two sublat- 
tice structure is only 0.47^s at T — > in zero field, 
as detected by elastic neutron scattering^. This quan- 
tum spin reduction indicates strong influence of quan- 
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turn fluctuations on the ground state. From the obser- 
vation of increase of ordered spin component in exter- 
nal field Tsyrulin et al^ conclude, that the fluctuations 
are suppressed by magnetic field. A related evidence 
of fluctuations suppression is the significant growth of 
Tjv in magnetic field, confirmed by neutron scattering 
and specific heat measurements^. A gap of Eq ~ 0.2 
meV in the spin-wave spectrum, detected by inelastic 
neutron scattering^, was ascribed to easy-plane (XY- 
type) anisotropy, keeping the spins within the be plane. 
The observation of minimum in susceptibility vs temper- 
ature dependence for a field directed perpendicular to be 
plane is consistent with quantum Monte-Carlo (QMC) 
simulations^, which predict the minimum of the suscep- 
tibility for HSLAF with a small XY anisotropy. 

We describe systematic investigations of 
Cu(pz)2 (0104)2 by means of multifrcquency electron 
spin resonance (ESR) spectroscopy and magnetization 
measurements for different orientations of magnetic 
field. Our main result is the observation and measure- 
ment of a weak in-plane anisotropy, not detected in 
previous measurements. This weak anisotropy induces 
remarkable features of the phase diagram. These are i) 
the spin-flop phase transition in a magnetic field applied 
along the easy axis, ii) a bicritical point and iii) a dip 
in Tn(H) dependence near the bicritical point. From 
antifcrromagnetic resonance spectrum we also find that 
the weak in-plane anisotropy is surprisingly changing its 
sign by a jump at the spin-flop point. Besides, this effect 
of abrupt anisotropy reversal arises as another phase 
transition at tilting the magnetic field, at a critical angle 
between the magnetic field and the easy axis. 



II. EXPERIMENT 

Samples of Cu(pz)2 (0104)2 have been grown in Clark 
University as described in Q. The lattice parame- 
ters of the monoclinic C2/c lattice are a = 14.072(5), 
b = 9.786(3) and c = 9.781(3) A; p = 96.458(4)°. Crys- 
tals are flat rectangular plaquettes colored blue, of the 
typical size of 2x2 mm 2 ; the plane of square lattice (6c 
plane) coincides with the plane of the plaquette. Sides 
of square crystal plaquettes are aligned at 45° to b and 
c axes, coinciding with the directing lines of magnetic 
square lattice, see a sketch on margins of FigfJ] 

ESR experiments were performed in Kapitza Institute, 
using a set of resonator spectrometric inserts in 4 He 
pumping cryostat with a cryomagnct. The frequency 
range from 5 to 140 GHz was covered. A spectromet- 
ric insert for 18 — 140 GHz range has a rotable sample 
holder, allowing to change the orientation of the sam- 
ple with respect to magnetic field during the experiment. 
A small amount of DPPH, free radical compound with 
g = 2.00, was used as a magnetic field label^. 

Magnetization experiments were performed at the De- 
partment of Low Temperature Physics and Superconduc- 
tivity of M. V. Lomonosov Moscow State University on 
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FIG. 2. (Color online) Left panel: resonance line temperature 
evolution in a sample, containing two types of domains, when 
i\) a = 15°. Right panel: resonance line temperature evolution 
in single domain sample when ip = 0° . In both cases £ = 0° , 
i. e. field lies in xy plane. The angles ip and £ are defined on 
Fig(4] A scaling factor of 0.25 is applied to 8 K line on the 
right panel. 
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FIG. 3. (Color online) Changes of resonance field and 
linewidth through Tat in a single domain sample for H || x at 
v = 27.46 GHz; corresponding set of ESR signals is shown at 
the right panel of Fig. [2] Lines are guide to the eye. 

Quantum Design 9 Tesla PPMS machine equipped with 
vibrating sample magnetometer (VSM) and at the Neu- 
tron Scattering and Magnetism Group in the Laboratory 
for Solid State Physics at ETH Zurich with the identical 
machine. The lowest available temperature was 1.8 K. 

III. ESR DATA 
A. Temperature evolution of ESR signal 

Magnetic resonance signal in Cu(pz)2 (0104)2 at tem- 
peratures above Tjy = 4.2 K corresponds to a typical 
exchange-narrowed paramagnetic resonance of Cu 2+ ions 
with anisotropic g-factor. The values of g-factor, ob- 
tained by high-temperature (T > 10 K) ESR measure- 
ments are g x = g y = 2.05 and g z = 2.28. A narrow 



3 



1.1 

1.0j 

H,0.9i 
T3 

Ll 

.2 

£ o.8| 
CO 0.9i 
1.0j 
1.1 



V=27.46 GHz B 90 

H_lz 



180- 




210 



middle a 



FIG. 4. (Color online) Polar plot of 27.46 GHz ESR field at 
rotating Cu(pz)2(C104)2 sample in the be plane; T — 1.3 K. 
Solid symbols correspond to actual experimental data, and 
open symbols are the repetition of solid symbols with the 
180° period. Symbols size corresponds to resonance intensity. 
Dashed lines are guide to an eye. A rough sketch of a two- 
domain sample and orientation of x,y,z axes within domains 
is shown. 
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FIG. 5. (Color online) Summary of angular dependencies 
in xy-plane for a domain. Solid lines are guide for the eye, 
dashed lines are theoretical calculation (biaxial model) with 
a shaded region around, marking a possible error due to the 
parameters uncertainty. Arrows indicate the field of param- 
agnetic resonance at corresponding frequency. 
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FIG. 6. (Color online). Evolution of the ESR line of a single 
domain sample at the rotation of the magnetic field within 
the xy-plane near the critical angle %p c — 10°, when regular 
mode 1/4 transforms to anomalous mode v a . Zero field cooling 
from 15 K to T = 1.3 K was performed before each record. 
Upper panel: v = 27.28 GHz. Lower panel: v = 37.7 GHz. 
Inset: Angular dependences of the intensities of ESR modes 
Vi (red circles) and v a (blue squares). Each half-open symbol 
corresponds to a line recorded after zero-field cooling. Solid 
symbols present ESR lines recorded without thermocycling, 
each rotation performed in a field 1.2 T. Thick lines are guide 
for an eye. 



Lorentzian line, with a halfwidth of about 5 ■ 10~ 3 T, 
broadens with cooling, and becomes unrcsolvable near 
Ttv- Below Tn the ESR response becomes strongly 
anisotropic. For H || z a broad signal transforms into 
a single narrow line, shifted from high-temperature po- 
sition, while for H || x,y two lines appear, as shown on 
left panel of Fig. [21 The resonance halfwidth shows a 
clear critical dependence near the phase transition tem- 
perature. The divergency in the line halfwidth together 
with the shift of the resonance position, as shown on Fig. 
[31 can be used as a marker of phase transition, allow- 
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FIG. 7. (Color online) Angular dependence of the field of 
49.31 GHz ESR in a two-domain sample of Cu(pz)2(C104)2 
at the temperature T = 1.3 K. Rotation is performed in the 
xz-plane of the domain A. Triangles — domain A, squares - 
domain B, circles - signal of the whole sample when domains 
become indistinguishable, crosses - T = 8 K ). Solid line 
presents the calculated low-temperature resonance field (bi- 
axial model)for i]) = 90°, dashed line - for rj) = and dotted 
line - theory for high-temperature paramagnetic resonance. 
Shaded region marks error boundaries of model calculation 
(see text). 



ing us to extract Tjv from the ESR data. The value of 
T/v=4.2±0.1 K is in agreement with the results of mag- 
netization measurements, as shown on the phase diagram 
(Fig. EH). 



B. Antiferromagnetic resonance 

The anisotropy within the be plane below Tjv results in 
the angular dependence of the resonance field, as shown 
on Figs. [U [5J This reveals two kinds of resonances with 
the identical rosette-like angular dependencies, which are 
shifted for 90° on Fig. |4j The relation between the in- 
tensities of these two kinds of signals is different for dif- 
ferent samples. This observation indicates a presence of 
two kinds of domains. The ratio of intensity of signals 
from two kinds of domains has the same value for zero- 
field cooling and field-cooling of the sample in the field of 
6 T, as well as at thcrmocycling through Tjv. Therefore 
we conclude, that these domains are crystallographic do- 
mains, for which b axes are rotated for 90°. We denote 
domains with orthogonal b (c) axes as domain A and 
domain B. One of the samples has the intensity for one 
rosette much stronger than for another one. For this, ap- 
proximately single domain sample, the relation between 
the volumes of domains of different types may be eval- 
uated, e.g., from the relation of intensities of ESR lines 
presented on FiglHl upper panel, between 0.82 and 0.88 
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FIG. 8. (Color online) Collection of ESR lines at T = 1.3 K 
for fields along tp = 0° (red solid) and tj) = 15° (blue dashed). 
Vertical dashed line denotes critical field H c = 0.42 T at this 
temperature. 



T. Such an estimation gives the number of spins, belong- 
ing to domain A, approximately 30 times greater than to 
domain B. This sample remained approximately a single 
domain one at numerous cycles of cooling from the room 
temperature. 

The rosettes shown on Fig|4] demonstrate a smooth 
evolution of the resonance field with the angle in the 
whole angle range except for the narrow range in the 
vicinity of 6-direction. This direction was identified for 
the nearly single domain sample by room temperature 
X-ray diffraction. At the angle -0 = 10° there is a step- 
like jump of the resonance field, shown in Figf?]and Fig. 
[5l Near the exact orientation of the external field along 
6-axis, i.e. when tilting ip does not exceed 10°, the reso- 
nance field is shifted to much lower field and this position 
can not be extrapolated from the smooth angular depen- 
dence in the main part of the field range. Therefore we 
denote the resonance observed at | ip \< 10° as anoma- 
lous mode v a . The redistribution of the intensity from 
the regular to the anomalous mode at a slow rotation 
of the field is shown on Fig. [6l One can see here that 
the transmission of the intensity between the two types 
of resonances has a character of a switching, it is per- 
formed within an interval of about 1° which may be a 
measure of the mosaic of the sample. Thus, a narrow 
phase transition at the angle variation is observed. 

We didn't observe any difference between the reso- 
nance signals, indicating this transition, when passing 
the critical angle in a field or in zero-field, and for pass- 
ing the critical angle at the temperature above or below 
Tjv- Tilting the field within the xz-plane conserves the 
anomalous mode, as shown on Fig. [71 at least in the range 
| £ |< 30°, where the difference between the anomalous 
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FIG. 9. (Color online) AFMR spectra at T = 1.3 K for 
field along easy axis. Solid line is the model calculation, 
dash-dotted line is the empirical formula |T]). Insert: shift 
from paramagnetic resonance frequency vs. field for direc- 
tions tp = 0° (empty circles) and tf> = 15° (solid squares) in 
xy plane. Solid and dashed lines are theoretical calculations, 
corresponding to these cases. Dash-dotted line is the same, 
as on main plot. 
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FIG. 10. (Color online) AFMR spectra for several field di- 
rections in xy plane and along z axis at T — 1.3 K. Solid 
lines are calculation according to (|B2[) with A y = 11 GHz 
and A z — 35 GHz, dashed lines are paramagnetic resonance 
with g xy — 2.05 and g z — 2.28. 



mode and an extrapolation for a regular mode may be 
detectable. The anomalous mode was observed only at 
/j,qH > 0.4 T. Below this field the resonance positions at 
-0 = 0° and ip = 15° are almost identical, as one can see 
on upper and lower records of Figj8] and on the low-field 
part of frequency-field dependencies on the inset of Fig. 
[9l At the same time, at hqH > 0.4 T there is a jump- 
like evolution of the resonance field and frequency in this 
range of angles. 

Further, we measured ESR fields for a set of frequen- 
cies (see examples of records on FiglHJ) at three principal 
directions of the magnetic field and at a tilting angle 
ip = 15°, as well as for two intermediate orientations in 
be plane. The corresponding frequency- field dependen- 
cies are presented on Figs. 191101 From these data we 
conclude, that the spectrum of frequencies of the anti- 
ferromagnetic resonance has two energy gaps, approxi- 
mately equal to 35 and 10 GHz and two branches in a 
magnetic field. For the direction of the magnetic field 
near 6-axis there is a mode softening at approaching the 
field of 0.42 T from the zero-field side. At fi H > 0.42 
T we observe the softened mode in the angular range of 
the regular mode and the anomalous mode in the narrow 
angle range | tjj |< 10°. By changing the angle i(j across 
the critical value toward -0 = 0, the ESR frequency is 
transposed from the value below the paramagnetic res- 
onance frequency g x [iB[i-QH jlixti to the value above it. 
This transposition is marked by an arrow on the inset of 
Fig. IH1 it occurs by a jump at crossing the critical angle 
rp c = 10°. This jump corresponds exactly to the jump of 
the resonance field, shown on Fig. [5J 



We note here that the observed frequency-field depen- 
dencies for field orientation in the whole solid angle, ex- 
cept for the range of the anomalous mode, may be well 
described by the calculated frequencies of the two sub- 
lattice antiferromagnet with a biaxial anisotropy and the 
easy axis directed along b, see, e.g. RefHol The calcu- 
lated frequencies are given in the Appendix A and pre- 
sented on Figs. [21 HD] by solid lines. Eight curves shown 
here, and the calculated angular dependencies shown on 
Fig. [7] are parameterized only by two energy gaps and 
three g-factors. The g- factors g x , g y , g z are measured 
independently in the paramagnetic phase and are not fit- 
ting parameters. Below the critical field of 0.42 T, the 
frequencies in the whole solid angle range of the magnetic 
field directions are described with that model. In partic- 
ular, a mode softening at H || b indicates the spin-flop 
transition. By this observation we can conclude that b 
is the easy axis direction. For the anomalous mode, ob- 
served at | V |< 10°, | £ |< 30°, fj, H > 0.42 T we use 
the empirical relation 

with A Q = 14 GHz at T = 1.3 K. This relation repre- 
sents the observed frequency at the unexpected position 
above (and not below) the paramagnetic resonance fre- 
quency at H > H c . 

Thus, the ESR data reveal a weak magnetic anisotropy 
in be plane and a spin-flop transition, as well as the 
anomalous mode v a appearing in the narrow angular 
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FIG. 11. (Color online) Isothermal low-field magnetization 
curves of Cu(pz)2(C104)2, field along x. An offset of 1.5 
emu/mol per curve is present. Upper insert: an expanded 
region around the temperature, where magnetization jump 
disappears. Lower insert: a few dM/dH curves are shown, 
allowing to locate phase transitions between disordered and 
ordered phases, marked by arrows. 



range of the field direction instead of a regular resonance 
of a biaxial antiferromagnet. 



IV. MAGNETIZATION 
A. Field along the easy axis 

The main feature of low-temperature magnetization 
curves at H || b is the presence of jump in magnetiza- 
tion corresponding to the spin-flop transition, detected 
by ESR. As it is shown at Fig. [TTJ magnitude of the 
jump increases with cooling, and its position shifts to 
lower fields. Jump in magnetization disappears around 
4±0.05 K, which is lower than T/y The sharp increase in 
magnetization can also be seen in M(T) curves, present 
on Fig. [TJ] Crossing the spin-flop phase boundary by 
temperature in a constant field also gives a very pro- 
nounced step in magnetic moment. This step disappears 
at /j,qH = 0.74 T, and above this field there appears a 
minimum on the M vs T dependence. Below the tem- 
perature of the minimum of the magnetization there is a 
kink, marking the onset of long range order. The min- 
imum and the kink are marked on Fig. 1121 Both min- 
imum and kink shift upwards in temperature with in- 
creasing field up to 9 T, though the former becomes less 
pronounced. Note, that there is no offset on Fig. [12J and 
stacking of the curves reflects the nonlincarity of magne- 
tization process. 

We derive the ordering point by a peak in the derivative 
d(MT) /dT, as suggested by Fisher—. We can also detect 
this transition by a peak in the derivative dM/dH of 
isothermal magnetization curve, as displayed at the insert 
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FIG. 12. (Color online) Normalized magnetization M(T)/H 
for various magnetic fields, directed along x. Insert: an ex- 
panded region around the point, where magnetization jump 
vanishes. 
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FIG. 13. (Color online) Phase diagram of Cu(pz)2(C104)2, 
field along easy axis. AF is 1 || x collinear antiferromag- 
netic phase, SF is spin-flop antiferromagnetic phase and PM 
is paramagnetic phase. Circles are features in M(T) curves, 
squares are features in M(H); red points correspond to tp = 0° 
orientation, and black — to ip = 12°. Lines are guide for an 
eye. In the insert an expanded region around the bicritical 
point is shown. 



of Fig. [TjJ] A final phase diagram, with points on phase 
boundaries obtained by both M(T) and M(H) scans, is 
present at Fig. UM 



B. Field in be plane 

On Fig. [TH a collection of magnetization derivatives 
dM(H)/dH for various directions of magnetic field in xy 
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FIG. 14. (Color online) Magnetization derivative dM/dH for 
various directions of magnetic field in xy plane at T — 2 K. 
Dashed lines show zero for corresponding curves. Inset: tem- 
perature dependencies of magnetization derivatives for direc- 
tions i> = 8° ± 2° and ip = 10° ± 2°. 
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FIG. 15. (Color online) Normalized magnetization M(T)/H 
for various magnetic fields, directed along z (left panel) and y 
(right panel) . Values of magnetic field between 0.3 and 1 Tesla 
are the same, that on Fig. 1121 The kink and the minimum on 
the curves are marked by a droplet and a cross, respectively. 



locate H c and T/v in the same way as was described be- 
fore for H || x. The resulting phase boundaries are also 
shown at Fig. [131 and the difference between phase di- 
agram for exact and misaligned orientations along x is 
observed only near the bicritical point. 



plane is present. dM(H)/dH curves show a sharp peak 
at ip — ► 0) which broadens with misalignment. Left wing 
of the peak is always smooth, in contrast to a discontinu- 
ity in the right wing of the peak. The discontinuity exists 
up to a critical angle ip c , observed in angular dependences 
of ESR, see Figs. IH [5] When tilt angle exceeds ip c , the 
dM/dH curve becomes completely smooth. The transi- 
tion between smooth and discontinuous types of deriva- 
tive, which affects mostly the right wing, occurs abruptly. 
This difference between the curves, corresponding to field 
tilts below and above ip c , is pronounced in a whole tem- 
perature range where spin-flop transition takes place, as 
shown in the inset of Fig. [T?] The discontinuity of differ- 
ential susceptibility dM/dH occurs exactly at the same 
magnitude and in the same angular range | ip \< ip c , as 
the anomalous ESR mode v a . 

With further increase in ip peak becomes less pro- 
nounced and almost disappears when ip approaches ~ 
45°. Magnetization curve at ~ 45° (this is the direc- 
tion along natural crystal facets) doesn't show a step, 
and demonstrates a smooth slope increase, as observed 
in Ref. El 

We have also performed a study of a phase diagram 
with magnetic field, slightly tilted from x axis. The 
orientation we chose was ip = 12 ± 2° > ip c . Here 
Cu(pz)2 (0104)2 still demonstrates increase in magneti- 
zation near H c , but the transition is regular, i.e. with 
a smooth derivative dM/dH on both sides of H c . We 



C. Field along middle and hard axes 

Curves of normalized magnetization for the fields, di- 
rected along middle axis y and hard axis z at Fig. [15] 
don't show spin-flop, but demonstrate an increase of Tn 
in a magnetic field. The only qualitative difference be- 
tween this two sets of curves is in the onset of minimum 
of M(T): for field along y the minimum appears for the 
applied fields between 1.5 and 2 T, while for field along 
z it is present even at H — > (as independently con- 
firmed by zero-field ac magnetization measurements^). 
Another feature of M(T) curves at H | y is the existence 
of inflection point below T/v at low fields, in contrast to 
the case of H || z. As the field is increased, the inflection 
disappears. 

V. DISCUSSION 

A. Biaxial model and spectra 

In section IIIII we have presented AFMR spectra for 
various field directions. These spectra follow biaxial 
collinear antiferromagnet paradigm for all magnetic field 
directions, except for a small solid angle corresponding 
to anomalous mode. The anomalous mode is observed in 
a solid angle of about of 10 _2 47r, close to easy axis, and 
only above H c . Resonance frequencies in this range of 
fields and angles correspond to a gapped branch (fjrj with 
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A a = 14 GHz. This unexpected effect can be described 
as an in-plane anisotropy switching caused by spin flop, 
i.e. the resonant frequencies are that of a two sublatticc 
biaxial antiferromagnet, for which x turns abruptly from 
the easy- into middle axis and y turns into easy axis at 
the spin-flop point. This conclusion is made on the base 
of the experimental observation of the anomalous mode, 
which appears at H > H c and has the frequency fol- 
lowing the relation (A3), corresponding to middle axis 
orientation of the field, instead of expected (A7) , derived 
for the easy axis orientation. 

We consider the magnetoelastic hypothesis, which 
might explain the switching of anisotropy at the spin flop 
point. The in-plane anisotropy, marking the easy axis, 
originates from rhombic distortion of square lattice, for 
Cu(pz)2 (0104)2 this distortion is due to a relative differ- 
ence of about ~ 10~ 4 between the lattice constants b and 
c. In principle, the antiferromagnetic ordering may cause 
a striction of the same order of magnitude^. Because 
of the magnetostriction, the in-plane anisotropy may be 
dependent on the magnitude and the direction of sub- 
lattice magnetizations, and, therefore, it should change 
when the spin flop transition takes place. One could ex- 
pect the y axis to become easy axis, and x to become 
middle axis immediately after the spin-flop. 

Nonetheless, a simple quantitative formulation of this 
approach, described in details in Appendix [B] does not 
capture a step-like angular dependence of AFMR field in 
xy-plane, and is in a contradiction with the observed re- 
lation between the zero-field gap and a critical field. The 
analysis of a complete Lagrangian, allowed by symme- 
try, should include several dozens of magnetoelastic and 
clastic terms and was not performed. 

Another possibility, presumably explaining the nature 
of anomalous mode v a , is the existence of a phase, other 
than collinear for H > H c , if; < ip c . This implies destabi- 
lizing of collinear phase by frustration when external field 
compensates in-plane anisotropy. However, our measure- 
ments do not support this hypothesis: magnetization 
curve M(H) and phase boundary Tn(H) are indistin- 
guishable for ip < tp c and ip > ijj c in fields above H c . 

The influence of a change of the direction and magni- 
tude of zero point fluctuations at the spin flop may be 
also of importance, because anisotropic spin fluctuations 
also contribute to the energy of anisotropy. 

Nevertheless, the nature of the anisotropy switching 
remains unclear. 

To give a connection with the previous work Q we 
derive a relation between zero-field gaps A VjZ of antifer- 
romagnetic resonance and the parameters of microscopic 
model Hamiltonian. In nearest-neighbor exchange ap- 
proximation the complete biaxial Hamiltonian reads as 

UA') i 

, ~ ~ ~ ~ x ( 2 ) 

-£ (5JyS?s?, + 5j g s;s;,), 
(i,i>) 



where 8 J y and S J z are parameters of so-called 'exchange 
anisotropy'. According to linear spin- wave approxima- 
tion, which have proven to be good for describing the 
fc-dependence of the spectrum in the vicinity of Brillouin 
zone center of Cu(pz) 2 (0104)2 @, the energy gaps are 
related to exchange anisotropy parameters as 

A y ,, = 2y/2J6J y , z . (3) 

In this spin- wave approximation the sublatticc magne- 
tization is supposed to be [Xb per magnetic ion, which 
is not the case of the Cu(pz) 2 (C104)2, where a strong 
quantum reduction of about 50 % is observed. A finer 
estimation for the case of S = 1/2 SLAFM, considering 
1/5 corrections, was given by Weihong et alM: 



A y>z ~ 1.2y/2JSJ y , z . (4) 

Thus, this equation may be used for an estimation of 
SJ VtZ . Spectroscopic gaps of 35 ± 2 and 11 ± 2 GHz are 
A 2 '= 1.68 ± 0.1 K and A y = 0.53 ± 0.1 K correspond- 
ingly. Hence, from (Q| we extract 5J Z = 53.2 mK and 
SJy = 5.3 mK. This corresponds to relative exchange 
anisotropy SJ Z /J = 3.1 • 10" 3 and SJ y /J = 3.1 • 10~ 4 . 
This is in agreement with previous neutron data, except 
for the parameter 5J y , which was not resolved by neutron 
scattering experiment. Wc can characterize the observed 
anisotropy switching, in terms of changing of parameters 
of Hamiltonian ((2]). It corresponds to transformation of 
SJ y into SJ*, which is of negative sign and equals —6.7 
mK. 



B. Phase diagrams 

Phase diagrams for x, y and z directions of magnetic 
field are present at Fig. 1161 For y and z directions 
the phase diagram are analogous, with monotonous in- 
crease of Tn- For field along x phase diagram is more 
complicated, with a bicritical point, where spin-flop, or- 
dered and paramagnetic phases meet. The phase dia- 
gram presents the spin flop transition and the bicritical 
point in addition to the phase boundaries reported in the 
previous work using neutron scattering and specific heat 
measurements^. 

The field, at which the antiferromagnetic resonance 
mode 1/3 is observed, also marks the spin-flop transi- 
tion (sec Appendix). From ESR experiment we get 
H Hf SR = 0.45 T at T = 1.3 K (sec Fig. ©, this field 
increases with temperature. Temperature dependence of 
H c derived from ESR is consistent with the magnetiza- 
tion measurements as shown on Fig. 1161 

Minima of M(T) are also plotted, showing different 
behavior for all three directions. While for H || z mini- 
mum persists up to H = limit, for different directions 
it appears only at some finite field, at which T m i n reaches 
T N . 
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FIG. 16. (Color online) Phase diagrams for H || x (upper 
panel) and H \\ yz (lower panel) directions. Points are exper- 
imental data: circles for M(T) and squares for M(H) features. 
Stars and diagonal crosses mark phase transitions determined 
by ESR, solid triangles are data from 0]. Horizontal crosses 
(x, y) and snowflakes (z) are minima in M(T). Solid lines are 
guide for the eye, dashed lines are fit ((5]l for T m i n (H) The fit 
for all three orientations is separately shown in the insert in 
the upper panel. 



The reason of this minimum may be qualitatively ex- 
plained by the following consideration: for an easy-pane 
antiferromagnet it is natural to have an anisotropic sus- 
ceptibility, which is larger for out-of-plane direction. In 
case of 2D antiferromagnet with 5 J z <§; J one should ex- 
pect this anisotropic behavior to rise only at low temper- 
atures, when T < J. Numerical simulations of HSLAF 
with a weak easy-plane anisotropy^ show, that tendency 
for increase of magnetization at H||z due to the onset 
of planar correlations overcomes the tendency for its de- 
crease due to short-range AF order. Thus a characteristic 
minimum in x(T) marks a crossover from Hcisenberg to 
XY behavior. Cuccolli et a££, using a QMC data analy- 
sis, suggested a formula for estimation of T m j„ in a case 
of easy-plane HSLAF model, 



4irp s 



hi 



(5) 



Here p s ~ 0.22 J is the renormalized spin stiffness, 5 e ff 
is the relative anisotropy and C ~ 160 is the dimension- 
less constant. It has also been found, that presence of 
external magnetic field in 2D magnets makes them effec- 
tively easy-plane and induces both Berezinsky-Kosterlitz- 
Thouless transition at finite temperature and a minimum 
in x{T) above it [TH]. In absence of long-range order, 
when a local order parameter 1 is formed, the orienta- 
tion 1 1 H provides an energy gain. Hence, the short- 
range order parameter becomes 2D instead of 3D in the 
isotropic case and the effective anisotropy energy in this 
field-induced XY behaviour is proportional to H 2 . 

For the orientation of the magnetic field H _L z for 
a strong enough field (gfiB/J-oH > \/6J z J), we consider 
'effective' easy-plane anisotropy induced by an external 
field. The easy plane of this anisotropy is perpendicular 
to the field. We take this anisotropy in the form derived 
in Rcf. M for HSLAFM 



k B J 



(0) 



where (3 is a dimensionless parameter. Here we 
disregard smaller anisotropy SJ y , as the experimental 
T m i n (H) in x and y directions is the same within the 
error bars. 

For the case H || z the easy-plane anisotropy originates 
due to the combination of the natural and field-induced 
anisotropy. We empirically combine these two factors 
which were analyzed separately in QMC simulation ;} 15 ' 16 



5Jz 
J 



k B J 



(7) 



Fitting experimental data for T m i n with equation ([5]). 
where S e ff is set as or S^ff, and parameters J — 

18.1 K and C = 160 are fixed, we yield p s ~ 0.24 J, SJ Z ~ 
0.023 K and (3 ~ 0.1, which is quite close to the result of 
Cucolli et al [l5j . Fits are shown on Fig. [THlwith dashed 
lines. This result can be considered as another indication 
of the 2D correlations developing in Cu(pz) 2 (0104)2 at 
T > Tjv. Nonetheless, value of SJ Z obtained by this fit is 
in better agreement with estimation by Eq. [3l which does 
not take into account quantum rcnormalization of the 
gap, than with Eq. |4l which considers 1/5 corrections. 

For H || x (i.e. field along the easy axis) there is a bi- 
critical point (H C ,T C ). Three phase transition lines meet 
in this point: second order paramagnetic to collinear an- 
tiferromagnetic phase transition (PM-AF), second order 
paramagnetic to flopped antifcrromagnctic phase transi- 
tion (PM-SF) and first order spin-flop phase transition 
AF-SF. In the vicinity of bicritical point the following 
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scaling equations are expected^: for AF-SF transition 
temperature dependence for critical field is 

H 2 (T)-H 2 c =A(^--lY (8) 

while for ordering transitions to AF and SF phases re- 
lations between magnetic field and ordering temperatures 
are 

H\T ) -Ht=A(L c - l y BAF (L- l y (9) 

and 

H\T)-HI = a{^--1^ +B SF (Jf--l) (10) 

correspondingly. For 'classical' 3D antifcrromagnet 
scaling exponent is known to be = 1.25 for uniaxial 
and <p = 1-175 for biaxial anisotropy. Theory also sug- 
gests amplitude ratio Q = Bsf/Baf = 1 for the former 
cas o 17 ' 18 . In contrast, for pure 2D case with easy-axis 
anisotropy bicritical point is expected^— to occur only 
at T — — a simple argument for that is the follow- 
ing: when easy-axis anisotropy is compensated by the 
external field, the system becomes equivalent to non- 
perturbed two-dimensional Hcisenbcrg model, which can 
possess long-range order only at zero temperature. The 
PM-AF and PM-SF phase boundaries, which meet at 
T = 0, are defined by 

|ff 2 (T)-H c 2 |cxT- 2 exp(-^. (11) 

The above equation is valid only in the absence of ad- 
ditional anisotropics and interlayer couplings, while in 
case of Cu(pz)2 (C104)2 both of these perturbations are 
present and the bicritical point is at T > 0. A numeri- 
cal proof for the latter statement can be found, e.g., in 
Monte-Carlo study of classical anisotropic XY antifcrro- 
magnet on square lattice [22[, where the phase diagram 
strongly resembles that of a 3D easy-axis AFM. Hence, 
phase diagram of Cu(pz)2 (0104)2 turns out to be an in- 
termediate case between ideal 2D and conventional 3D 
anisotropic antifcrromagncts. Straightforward fit of ex- 
perimental data with equations (|8|9|10|) gives bicritical 
point at T c = 3.97 K, fioH c = 0.73 T with scaling ex- 
ponent = 1.4 and amplitude ratio Q = 1.78. Fixing 
the value of <j> to theoretically suggested for 3D antifcrro- 
magnet fort = 1.175 leads to T c = 3.99 K, fioH c = 0.738 
T and Q = 1.54, but with a worse fit quality. Data 
and fits are present at Fig. 1171 together with numerical 
quality criterion — sum of average least squares for all 
three formulas ()8I9I10[) . It can be concluded that reli- 
able estimation of universal parameters from our data is 
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FIG. 17. (Color online) Scaling in the vicinity of bicritical 
point. Blue solid and red dashed lines are fits (|8I9I10[) with 
cj> — 1.4 and <f> = 1.175 correspondingly, dotted lines are ex- 
tension of JS} to the T > T c region. Left insert: magnified 
part of main plot. Right insert: fit deviation (see text) versus 
critical exponent <j>. 



cf> = 1.4 ± 0.2 and Q = 1.8 ± 0.2, and bicritical point is 
located at (j, H c = 0.730 ± 0.006 T and T c = 3.97 ± 0.03 
K. Region AT ~ 0.2 K where the scaling equations are 
fulfilled is about 5% of T c , which is significantly larger 
than for classical three-dimensional uniaxial antifcrro- 
magnet MnF 2 (AT/T C ~ 10~ 3 [3), though smaller than 
for quasi-2D compound Rb 2 MnF 4 (AT/T C ~ 0.26 [H|), 
with a purely uniaxial anisotropy. These facts, as well 
the larger value of critical index <f> = 1.4 > fop result in 
a conclusion, that Cu(pz) 2 (0104)2 presents an intermedi- 
ate behaviour between 3D and 2D models in the vicinity 
of bicritical point. 

The observed dependence of the phase diagram in the 
bicritical point range on the field orientation is nat- 
ural, because bicritical point is very sensitive to field 
misalignment^. On the contrary, away from the region 
around H c , the phase boundaries for tp = 0° and if) = 12° 
coincide within error bars. 



VI. CONCLUSIONS 

In the present work we have studied AFMR spectra 
and magnetization curves of HSLAFM Cu(pz)2(C104)2- 
These measurements reveal the presence of biaxial 
anisotropy in Cu(pz)2 (0104)2, instead of easy-plane for- 
mulation used earlier. From the ESR experiments we 
have derived two energy gaps A z ~ 35 and A y ~ 11 
GHz. The weak in-plane anisotropy is responsible for 
the spin-flop phase transition at \iqH c — 0.4 T in H || x 
direction. The AFMR spectra also show, that weak in- 
plane anisotropy is changing its sign at the spin-flop tran- 
sition. This anisotropy reversal, occuring in a manner of 
switching, may be also observed as a phase transition at 
changing the orientation of the magnetic field within the 
be plane, at the critical angle 10° with respect to the easy 
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axis direction. The conjecture that this anomaly might 
be of magnetoelastic origin may, by simplified treatment, 
explain the anisotropy reversal at the spin flop, but is not 
consistent with the step-like angular dependence of ESR 
field or frequency. It is also not consistent with the ob- 
served relation between the energy gap A y and spin-flop 
field H c . The nature of the abrupt reversal of the weak 
anisotropy remains unclear. 

The hypothesis which is probably worth to analyze the- 
oretically is a possible change of the direction and mag- 
nitude of zero point spin fluctuations at the spin flop. A 
change of the contribution of fluctuations to the energy of 
the anisotropy may also change the effective anisotropy 
of the ordered spin component. 

The field-dependence of the temperature of the mini- 
mum on M(T) curves for three principal orientations are 
found to be in agreement with the results of numerical 
simulation of HSLAFM^. The increase of Tn in external 
field has been found for all orientations in agreement with 
previous measurements. Scaling exponent <f> = 1.4 ± 0.2 
of phase boundaries near the bicritical point is interme- 
diate between 2D and 3D models. 

Accounting for the observed weak anisotropy might be 
significant for correct estimation of other weak interac- 
tions, e.g. next-nearest neighbor and intcrlayer exchange 
from experimental data3i Similar anisotropy can be 
present in another HSLAFM's of Cu-pz family (namely, 
Cu(pz) 2 (BF 4 ) 2 and [Cu(pz) 2 (N0 3 )](PF 6 )), as according 
to Xiao's magnetization data^ there are signatures of 
spin-flop transitions as well. 

We note that for deuterated Cu(pz) 2 (C104)2 the dif- 
ference between b and c, resulting in the weak anisotropy, 
is larger than for a regular sample with hydrogen, so it 
would be of interest to test if the in-plane anisotropy is 
the same in a deuterated sample. Other possible future 
experiments include NMR and neutron scattering with 
field along x to probe the magnetic structure as well as 
search for anomaly in magnetic field dependencies of low- 
temperature elastic properties (magnetostriction, ultra- 
sound propagation etc). 
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Appendix A: AFMR frequencies of a biaxial 
antiferromagnet 

A theory for AFMR in a two-sublattice antiferro- 
magnet with biaxial anisotropy has been developed in 
1950's^. For both orientations of the magnetic field 
along hard- and middle-axis (z and y), in the ground 
state we have the antifcrromagnctic order parameter 
1 j| x, and this orientation of 1 is independent of external 
field magnitude. 

Magnetic resonance frequencies for H || z are 

v zz = ^(|gVff) 2 + Ai, (Al) 
v zy = A y ; (A2) 
and for H j| y we have 

v yy = y(f^W) 2 +A2, (A3) 
v v , = A z . (A4) 

At the orientation of the magnetic field along the easy 
axis the case is more complicated, as the ground state is 
field-dependent. There is a spin-flop transition with an 
abrupt change from 1 || x to 1 || y. The critical field of 
this transition is 

Mo# c = 2tt-^ (A5) 

9x^B 

This transition is accompanied by a jump in magnetiza- 
tion. The ESR frequencies below and above H c are the 
following: 
For H < H c : 



(A6) 

I 

The resonant mode v% with the vertical vz(H) depen- 
dence corresponds to spin-flop transition, as the system 
is allowed to absorb energy in a band of frequencies in the 
critical point. Modes v\ and v± are softened at the crit- 
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ical field H c . For the intermediate field orientations we 
calculated the frequencies of spin resonance numerically 
within the same formalism. 



Appendix B: Magnetoelastic correction 

For description of ESR modes at T — > we use macro- 
scopic exchange symmetry formalism^. This formal- 
ism, in particular, reproduces the results of a mean-field 
theory of a two-sublattice antiferromagnet with biaxial 
anisotropyi^. In the framework of the exchange ap- 
proach the spin structure is considered to be collincar, 
and the anisotropy of a relativistic origin, and magnetiza- 
tion, induced by the external field, are taken as perturba- 
tions. Though being applicable only in fields H <C H 8at , 
this formalism is model-independent and allows easy in- 
troduction of additional anisotropy terms. As the satu- 
ration field in Cu(pz)2 (0104)2 constitutes almost 50 T, 
restriction on the field magnitude is not an issue for the 
exchange symmetry formalism applicability. Our calcu- 
lations are based on the following Lagrange function per 
mole of the compound (in COS system): 



XX 
2 7 2 



i + 7 [Hx 1] 



U a . 



(Bl) 



Here 7 



9(J>b 



is the gyromagnetic ratio, unit vector 



1 with the orientation, given by angles (p and 8 shown 
on Figf5J is the order parameter, H is magnetic field; 
X± is the magnetic susceptibility in the direction, per- 
pendicular to 1. Equation (|B1[) also implies x\\ = at 
zero temperature. Term U a = f]l y + Ql is the anisotropy 
energy. We assume positive constants r\ and C, V < C 
Hence 1 || x minimizes anisotropy energy. 

The ground state and magnetic resonance frequencies 
may be calculated using this Lagrange function, as de- 
scribed in Refl24l The ground state and the spectrum 
are identical to that of Ref. [Io|, described above. The 
anisotropy constants may be expressed via energy gaps: 
7] = X ±(2^A 2 )/2 7 2 and C - x±(2^A 2 ) z /2 7 2 , where 
A z > A y are the energy gaps, which one actually ob- 
serves in the ESR experiment. With this substitution, 
Lagrange function (|B1|) is 



£ = ^Q(l + 7 [Hxl] N2 



(27tA z ) 2 )2 (2nA y )\ 2 
2 2 2 y 



(B2) 



and corresponding potential energy in non-zero mag- 
netic field is 



(2nA^ (2nAyl 
2 z+ 2 



(B3) 
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Magnetic field (T) 

FIG. 18. (Color online) Spectra of model (|B4|) for vari- 
ous values of quartic term. Field is directed along x. Ar- 
rows mark frequencies at which angular dependencies for 
( 2nhA v , 

are shown m the insertion; dashed lines 
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are the expectation for paramagnetic resonance. 
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FIG. 19. (Color online) Spectra of the model with a quartic 
term —j3lt, H || x. Dashed line is paramagnetic resonance. 



Monoclinic symmetry allows for another second-order 
term, l y l z . Such a term results in a tilt of hard and middle 
anisotropy axes, leaving easy axis undisturbed. In our 
experimental data, related to xz plane rotation of the 
magnetic field (Fig. [7]), we do not notice any significant 
tilt of middle axis from c direction, and, therefore, we do 
not take l y l z term into account. 

We have to note, that anisotropic term rjl 2 originates 
from a weak orthorhombic distortion of a square lattice. 
Due to this distortion, the lattice constants b and c dif- 
fer in a relative sense for about of (~ 10~ 4 ). Hence, 
this term should be small in comparison with, e.g. £Z 2 , 
and can be comparable with the contributions of higher 
order in components of 1. There is a term B(l x l y ) 2 
among the fourth-order terms, allowed by symmetry for 
Cu(pz)2 (0104)2- This term couples the components l x , l y 
and could result in the 'anisotropy reversal' as a result of 
1 reorientation. Indeed, considering a modified Lagrange 
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u=c/2 r + n/2 i v 2 + pi x 4 
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FIG. 20. (Color online) Spectra of the model with a quartic 
term /3lt, H || x. Dashed line is paramagnetic resonance. 
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FIG. 21. (Color online) Spectra of the model with a quartic 
term —/3ly, H [| x. Dashed line is paramagnetic resonance. 
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FIG. 22. (Color online) Spectra of the model with a quartic 
term ply, H || x. Dashed line is paramagnetic resonance. 



function 

C =^(l (i + 7[Hxl] N2 



(B4) 



we obtain approximate frequencies, corresponding to in- 
planc fluctuations of 1, for ground states before and after 
reorientation: 



iff"' 1 



and 



Hi II v) = 



The field-independent constants under the square root 
signs in this relations may be treated by use of model 
(B2) and relations (A6,A7) as anisotropy constant for 
in-plane anisotropy (note that these constants should be 
taken with the opposite signs). Thus, if B > is large 
enough, the in-plane anisotropy effectively changes its 
sign. Significant value of B might be provided by a mag- 
netoelastic interaction. There are 13 elastic terms of the 
form of UikU mn and 20 magnetoclastic terms of the form 
of u.ikl m l n allowed by symmetry for Cu(pz)2 (0104)2 [25j . 
From these terms we take for the magnetoelastic contri- 
bution to potential energy (|F33|) the following terms: 



and 



E. 



oc 



2 



(B5) 



(B6) 



These terms couple l x and l y . Here Uik are components 
of the strain tensor. 

Minimization of energy, including described above 
magnetoelastic correction (|B5IB6j) . with respect to strain 
variable u xy will result in the magnetoelastic correction 
in the form of Bill?,. 

We numerically calculate ground state and spectrum 
of the model (IB4I). For our numerical work we choose 



parameters A z = 37 and A y = 15 GHz, and follow 
the perturbations, introduced by quartic term Bl x l y . At 
Fig. [18] numerically calculated spectra for field along 

x and different values of B (namely, — ( ^ x j 1 ^ 

1 X ± \2irHA y/ 

0,0.1,0.5,1,1.5,2) are present. Indeed, we can choose 
the value of B, which is large enough to push mode 1/4 
above paramagnetic resonance, as it is clearly seen on 
Fig. 1181 But there remains a crucial difference between 
the properties of model (|B4|) and experimental data, be- 
cause the angular dependence is not described. We found 
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from ESR experiment that angular dependence of reso- 
nance line shows a step when switching from v a to 1/4, 
the anomalous mode v a exists in a narrow angle range 
near the x-axis, and outside of this range all the data is 
perfectly described by a simple biaxial model. In con- 
trast, the model with strong B term (|B4j) shows smooth 
angular dependence (see insert on Fig. I18|) . which differs 
significantly from both simple biaxial model (|B2[) and 
experimental data (as on Fig. [5]). Furthermore, as it is 
seen from Fig. [TSl enhancement of B term also increases 
a lower zero-field gap, while critical field of spin-flop tran- 
sition is still determined by Ay alone, as l^ly combination 
gives the same contribution to energies of both 1 || x and 
1 || y phases. Therefore, the magnetoelastic approach 
(|B4j) predicts the transposition of the antiferromagnetic 
resonance frequency above the value of the paramagnetic 



resonance frequency, which is one of the manifestations of 
the reversal of the in-plane anisotropy. Nevertheless, at 
the same time, the angular dependence of the resonance 
field and relation between the gap and critical field do 
not correspond to the experiment even qualitatively. 

Other related quartic terms also could not describe the 
observed anomaly. We have analyzed in the same way the 
influence of anisotropic terms ±/3Z^ and ±/3^. The re- 
sults of the calculation of the ground states (equilibrium 
values of ip, 9) and frequency-field dependencies arc given 
in Figs. Uni HOI HH I2U One can see, that it is impossible 
to find a value of (3 which would correspond to the ob- 
served pulling of the frequency above the paramagnetic 
resonance frequency for H > H c along with the soften- 
ing of the mode at H < H Cl and with the valid relation 
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